Unfolding first-principles band structures 
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A general method is presented to unfold band structures of first-principles super-cell calculations 
with proper spectral weight, allowing easier visualization of the electronic structure and the degree 
of broken translational symmetry. The resulting unfolded band structures contain additional rich 
information from the Kohn-Sham orbitals, and absorb the structure factor that makes them ideal for 
a direct comparison with angular resolved photoemission spectroscopy experiments. With negligible 
computational expense via the use of Wannier functions, this simple method has great practical value 
in the studies of a wide range of materials containing impurities, vacancies, lattice distortions, or 
spontaneous long-range orders. 
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The electronic band structure is no doubt one of the 
most widely applied analysis tools in the first-principles 
electronic structure calculations of crystals, especially 
within the Kohn-Sham framework [1] of density func- 
tional theory [2]. It contains the basic ingredients to 
almost all the textbook descriptions of crystal properties 
(e.g. transport, optical and magnetic properties, and the 
semiclassical treatment [3]). Furthermore, the theoret- 
ical band structure, when formulated within the quasi- 
particle picture of the one-particle Green function, has 
a direct experimental connection with angular-resolved 
photoemission spectroscopy (ARPES). 

However, the usefulness of the band structure, as 
well as the agreement with ARPES spectra, diminishes 
rapidly when a large "super cell" is involved. The use 
of super cells is a common practice in modern first- 
principles studies when the original periodicity of the sys- 
tem is modified via the introduction of "external" influ- 
ences from impurities or lattice distortions. They are also 
widely applied in the presence of spontaneous transla- 
tional symmetry breaking, say by a charge density wave, 
a spin density wave, or an orbital ordering. As illustrated 
in Fig. [l] when the period of the super cell grows longer, 
the corresponding first Brillouin zone of the super cell 
(SBZ) shrinks its size. In turn, bands in the first Bril- 
louin zone of the normal cell (NBZ) get "folded" into the 
SBZ. For a very large super cell, the resulting SBZ can 
be tiny in size but contain a large number of "horizontal" 
looking bands that no longer resemble the original band 
structure or the experimental ARPES spectra, and cease 
to be informative besides giving a rough visualization of 
the density of states (DOS). The information is now hid- 
den in the Kohn-Sham orbitals, instead of the dispersion 
of the bands. 

In this Letter, by explictly utilizing these Kohn-Sham 
orbitals, we present a method to unfold the band struc- 
ture of the SBZ back to the larger NBZ with proper 
spectral weight. Making use of the corresponding Wan- 
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FIG. 1: (color online) Illustration of band folding in the super 
cell calculations: (a) band structure of a 2D one-band first- 
neighbor tight-binding model, (b) the same obtained from a 
4x4 super-cell calculation, and (c) the same obtained from a 
16x16 super-cell calculation. Pannel (d) shows the DOS. 

nier functions, the method can be greatly simplified to 
negligible computational cost. The resulting unfolded 
band structure incorporates explicitly the structure fac- 
tor and thus facilitates significantly a direct comparison 
with ARPES experiments. Furthermore, the unfolded 
band structure illustrates very clearly the influence of the 
symmetry breaker (e.g.: impurities, vacancies, dopants, 
lattice distortions) via direct comparison with the nom- 
inal normal-cell band structure. In the case of sponta- 
neous symmetry breaking, it gives a direct visualization 
of the strength of each band's coupling to the order pa- 
rameters. In light of the amazingly rich information, we 
expect countless applications of this simple method to 
a wide range of studies employing super cells, including 
systems with charge density wave, spin density wave, or 
orbital ordering, and in the studies of impurities and lat- 
tice distortions, to name a few. 

Theoretically, the folding of the bands results from the 
introduction of additional coupling, Vkj^kiji, between the 
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originally uncoupled Kohn-Sham orbitals \kj) and \krjr) 
in the NBZ. (Here k and j denote the crystal momentum 
and the band index.) This coupling extends the period 
of Kohn-Sham orbitals to a longer one compatible with 
the size of the super cell. Equivelantly, this coupling, 
no matter how small it is, mixes the original orbitals of 
different normal-cell crystal momentum k and forces us 
to label them with a supercell crystal momentum K as 
the new quantum number in the SBZ. (In the following, 
upper-/lower-case symbols refer to variables correspond- 
ing to the super/normal cell, respectively.) Our method 
is based on the simple idea that unless V is extremely 
strong, it is much more convenient and informative to 
represent the band structure or more precisely the spec- 
tral function A = — ImG/vr of the retarded one-particle 
Green function, G, not in the new eigen-orbital \KJ) ba- 
sis, but in the \kj) basis of the normal cell instead: 

where Gq represents a conceptual system with the pe- 
riod of the normal cell before V is applied. Clearly, G 
smoothly recovers the original period of Go as V ap- 
proaches zero. Thus, 

Ak,^kM = Y^\{kj\KJ)\''AKj,Kj{u^) (2) 

KJ 

should resemble the band structure of the normal cell 
with deviations in both the dispersion and in the spec- 
tral weight that reflect the effects of V . Note that 
while the coupling V introduces non-diagonal elements of 
■^kj,kiii{'^), we focus only on the diagonal elements here 
for simplicity, without loss of generality. It is straightfor- 
ward to show that in the case of F = 0, the weight of the 
bands that follow the bands of the normal cell is exactly 
one, and that of the rest of the folded bands vanishes. 
One thus recovers exactly the original band structure of 
the normal cell as expected. That is, the unfolded hand 
structure is invariant against any arbitrary choice of su- 
per cell. 

In addition, it is often desirable to also measure in 
each band the contribution of local orbitals with well- 
defined characters (e.g.: Px, dxz, Cg, or bonding/anti- 
bonding) . This can be achieved rigorously via the use of 
local Wannier orbitals \rn): 

Akn,kn{^) = ^\{kn\KJ)\''AKJ,Kj{^), (3) 
KJ 

where \kn) = \rn){rn\kn) = \rn)e^''^ /y/l are the 
Fourier transform of the Wannier orbitals \rn) of orbital 
index n and associated with the lattice vector r. (Here 
I denotes the number of k-points in the NBZ.) Given 
a consistent definition of the Wannier functions of the 
super-cell calculation that maps \RN) of the super cell to 
\R-\-r, nf) of the normal cell, where r = r{N) is a normal- 
cell lattice vector within the first super cell, and n/ — 



ni{N) is the corresponding normal-cell orbital index, the 
use of Wannier function also reduces dramatically the 
computational expense by turning the factor 

{kn\KJ) = ^{kn\RN){RN\KN){KN\KJ) 

RN 

= + r{N), n/{N)) {RN\KN){KN\KJ) 

RN 

RN 

= v^^e-*'^(^)5„,„,(^.)5[fe],K(i^iV|A'J) (4) 

N 

into merely a structure factor that is a sum of coeffi- 
cients of the eigen-orbital \KJ) of the super cell in the 
Wannier function basis, modulated by the proper phase 
that encapsulates the internal position in the super cell. 
Here [k] denotes the k-point folded into the SBZ from k. 
Since Akj,kj is just a delta function at the eigenvalue 
(5(a; — tKj), this final expression in essence requires only 
a simple coding to plot all the eigenvalues of the super 
cell in the larger NBZ with a proper weight. 

Of course, the above definition only makes sense when 
the Wannier functions \RN) ^ \rn) and \RNf) o |r/n) 
(that are translational symmetric in the normal cell unit: 
same n different r) are approximately identical. There- 
fore, the "gauge" [1] of constructing \RN) and \RNf) 
with the same n must be controlled accordingly. In the 
presence of a potential that breaks the translational sym- 
metry of the normal cell, for example, coming from a 
CDW, lattice distortions, impurities, etc., the commonly 
employed [SHZ] maximally localized Wannier function |4] 
and other minimization-based methods (SUHj risk defining 
the gauge differently in the super cell in favor of better lo- 
calization, and thus should be used with extreme caution. 
We found that a maximum projection method [T0HT2] 
with consistent projection between the normal and the 
super cells works well to satisfy this requirement. Equa- 
tions ([3| and Q should in principle also be applicable 
in many existing codes employing atomic center local or- 
bitals as basis [13 |T3], as long as the non-orthogonal 
nature of those bases is taken into account. Of course, 
these methods do not benifit from the energy resolution 
of the Wannier functions that allows unfolding only the 
bands within the physically relevant energy range. 

The unfolded band structure also has an important di- 
rect connection to the ARPES measurement. For sys- 
tems with enlarged unit cells due to weak symmetry 
breaking, the ARPES spectra typically shows different 
band structures in different Brillouin zones of the super 
cell, distinctly different from the results of first-principles 
calculations, which have all the bands in the SBZ. In 
some cases, the observed ARPES spectra might even ap- 
pear ignorant about the SBZ [TS]. This significant mis- 
match is typically regarded as the effect of the "matrix 
element" and left unaddressed by both theorists and ex- 
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perimentalists, making a direct comparison very difficult. 
Witfiin the "sudden approximation" , the ARPES inten- 
sity is proportional to [lH [17] 



J2\^-{f\p\KJ)\'AKJ,KjH 
KJ 

^ \e-{f\p\kn)\^\{kn\KJ)\^AKj,Kj{io) 

KJkn 

= Y,\e-{f\p\kn)\^AknM^), (5) 

kn 

where e denotes the polarization vector of light, and 
I/) the "final state" of the photoelectron. Clearly, ex- 
cept the polarization dependent dipole matrix element, 
|e- (/|p|fcn)p, the unfolded spectral function, Akn,kn{^)i 
contains almost the full information of the experimen- 
tal spectrum by absorbing the additional structure fac- 
tor {kn\KJ), absent in the typical super-cell solution, 
■^Kj,Kj{^)- Obviously, the inclusion of this additional 
matrix element would facilitate significantly the compar- 
ison between the theory and the ARPES experiment. 

As an example, let's consider the effect of Na impu- 
rities in Na-doped cobaltates, Naa;Co02 at a; = 1/3. In 
typical first-principles studies 119) . the impurity is in- 
corporated via a super cell as demonstrated in Fig. ^h) 
in comparison with the undoped normal cell shown in 
Fig. [2]^ a). Fig. [2|d) and (c) show the corresponding 
band structures obtained with standard DFT calcula- 
tions. Since in this example the super cell is three times 
larger than the normal cell, the corresponding SBZ is 
three times smaller and contains three times more bands. 
Clearly, even for such a small super cell, the change of 
the size/orientation of the SBZ and more importantly the 
large number of folded bands, make it practically impos- 
sible to cleanly compare with the band structure in the 
NBZ of the undoped parent compound. In fact, to many 
untrained eyes, these two band structures may appear 
entirely unrelated. 

By contrast, the unfolded band structure shown in 
Fig. ^e) , demonstrates a strong resemblance to the band 
structure of the undoped compound. This allows a clear 
visualization of the effects of the (periodic) Na impurities 
on the original Co and O bands. Specifically, besides the 
introduction of additional Na-s bands, one observes shifts 
in band energies, gap openings and the nearby "shadow 
bands" , all of which refiects the infiuence of the Na impu- 
rity on these bands. What is really nice here is the clean- 
ness of the unfolded band structure in general, owing to 
the weak intensity of the shadow bands. As expected, the 
influence of the Na impurity is only minor on most Co- 
d and 0-p bands, while the Na-s bands themselves show 
sizable effects of broken translational symmetry. The size 
of the gap opening and the intensity of the shadow bands 
actually reflect directly the strength of each band's cou- 
pling to the broken translational symmetry of the nor- 
mal cell (in this specific case, to the charge-density-wave 
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FIG. 2: (color online). Lattice structures of (a) C02O4 (nor- 
mal cell) and (b) Na2Co60i2 (super cell), the corresponding 
band structure of (c) the normal cell and (d) super cell calcu- 
lation, and (e) the unfolded band structure of the super cell. 
Insect illustrates the effects of weak translational symmetry 
breaking via spectral functions over the region [-4.6eV,-4.2eV] 
and [frAf, |rA/]. 



order parameter introduced by the periodic presence of 
Na atoms.) Of course, for a simulation of randomly po- 
sitioned impurities, these CDW-related features are en- 
tirely artificial, and the unfolded band structure makes 
apparent the alarming limitation of common practice of 
using small super cells in the study of impurities. On the 
other hand, in many other cases, for example the super 
modulation of the lattice, these features would actually 
correspond to a physical order parameter and provide 
valuable information. 

As another example, let's consider a spontaneous or- 
bital ordering in A- type anti- ferromagnetic LaMnOa. 
Figure [s] (c) and (d) show the similar comparison of 
band structures without and with the long-range stag- 
gered orbital order, corresponding to unit cells shown in 
Fig. [3] (a) (normal cell) and Fig. |3] (b) (super cell), re- 
spectively. Both results are obtained via the LSDA+U 
(C/=8eV, J/f =0.88eV) approximation without lattice re- 
laxation for simplicity, without loss of generality. By 
comparing the band structures with and without the or- 
bital order on the equal footing, the detailed information 
of the spontaneous orbital order should be visualized ex- 
plicitly. 

Just like in the Na2:Co02 case, the straightforward re- 
sults of the orbital ordered (00) band structure (Fig. [3] 
(d)) of the super cell calculation [5D] can hardly be com- 
pared with the non-00 one (Fig. [s] (c)). By contrast. 
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(a) Normalcell (La2Mn20e) (b) Supercell (La4Mn40i2) 
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FIG. 3: (color online). Lattice structures of (a) normal cell 
and (b) doubled-size supercell of A-type anti-ferromagnetic 
ordered LaMnOa, without and with the orbital-ordering, re- 
spectively, the corresponding band structure of the normal 
cell (c) and the super cell (d), and the unfolded band struc- 
ture of the super cell (e) indicating the orbital ordering gap, 
Aoo. Red and blue bands denote the and — orbital 
characters of the spin- majority channel, and green bands gives 
both Eg characters of the spin-minority channel. 



the unfolded band structure (Fig. [s] (e)) of the 00 case 
resembles strongly the non-00 case. In fact, one finds 
that only those bands of Mn-e^ character (red, blue, and 
green) show strong coupling to the 00 order parame- 
ter with large gap openings and intensive shadow bands, 
while the rest of the bands are basically uncoupled to the 
orbital order. In addition, from the significant energy 
gain corresponding to the large 00 gap (Aoo) near the 
Fermi level of the red and blue bands, it is apparent that 
essentially the orbital order is driven only by the spin- 
majority Eg-orbitals (red and blue). All these effects are 
of course entirely consistent with the existing "electronic 
interaction assisted Jahn- Teller picture" , in which 

the degenerate Mn-e^ orbitals split to gain energy and 
stabilize the system at low-temperature. However, this 
unfolded band structure represents probably the best vi- 
sualization of such physics in real materials with details 
of first-principles calculations. 

In conclusion, a simple method for unfolding first- 
principles band structures of super cell calculations is 
presented. Proper spectral weights are obtained with 
negligible computational cost by making use of the Kohn- 
Sham orbitals with the help of carefully chosen Wannier 
functions. The inclusion of the structure factor in the re- 
sulting unfolded band structure makes it ideal for direct 



comparison with the ARPES measurement. The result- 
ing unfolded band structures allow an easy visualization 
of each band's coupling to the order parameter of spon- 
taneous broken translational symmetry, as well as their 
couplings to the external symmetry breakers like the im- 
purities and lattice distortions. Our method should prove 
valuable in the study of a wide range of problems requir- 
ing the use of super cells, including systems with impuri- 
ties, vacancies, and lattice distortions, and broken sym- 
metry phases of strongly correlated materials, to name a 
few. 
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